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We study the transverse momentum dependence of the 
Landau-Pomeranchuk-Migdal effect in QED, starting from 
the high energy expansion of the solution of the Dirac equa- 
tion in the presence of an external field. The angular in- 
tegrated energy loss formula differs from an earlier expres- 
sion of Zakharov by taking finite kinematical boundaries into 
account. In an expansion in powers of the opacity of the 
medium, we derive explicit expressions for the radiation cross 
section associated with N = 1, 2 and 3 scatterings. We verify 
the Bethe-Heitler and the factorization limit, and we calcu- 
late corrections to the factorization limit proportional to the 
square of the target size. A closed form expression valid to 
arbitrary orders in the opacity is derived in the dipole approx- 
imation. The resulting radiation spectrum is non-analytic in 
the coupling constant which is traced back to the transverse 
momentum broadening of a hard parton undergoing multiple 
small angle Moliere scattering. In extending the results to 
QCD, we test a previously used dipole prescription by com- 
paring to direct pQCD results for A'^ = 1 and 2. For N — 1, 
the QCD dipole prescription reproduces exactly the Bertsch- 
Gunion radiation spectrum. For TV = 2, we find a sizeable 
correction which reduces to a multiplicative factor 17/8 at 
large separation. 



I. INTRODUCTION 

In QED, the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect interpolates between the Bethe-Heitler and factoriza- 
tion limit for the radiation spectrum of a charged particle 
undergoing multiple, N say, scatterings in a medium. If 
the separation between scattering centers is very large, 
then the radiation off these centers reduces to a sum of 
radiation spectra for N small-angle scatterings - this is 
the Bethe-Heitler limit. In the opposite limit, when the 
scattering centers sit too close together to be resolved 
by the emitted photon, the observed radiation factorizes 
into a product of a single scattering Bethe-Heitler spec- 
trum for momentum transfer q = J2i=i 1i ^^^ ^^^ elastic 
cross section for the momentum transfer q accumulated 
over TV small-angle scatterings. 

As first noted by Landau and Pomeranchuk [ |l|, |[ , the 
relevant length scale for the interpolation between Bethe- 
Heitler and factorization limit is the coherence length 
(formation length) If, determined by the longitudinal 
momentum transfer qj. 



h = yqi 



(1.1) 



This characterizes the longitudinal scale on which the 
radiated particle becomes distinguishable from its radi- 
ating parent. Scattering amplitudes for the radiation off 
different scattering centers interfer destructively if their 
separation is less than If : the coherent factorization limit 
is suppressed in the ultrarelativistic limit with respect to 
the incoherent Bethe-Heitler limit. 

For a quantitative description of the LPM interference 
effect, the relative phases of the different contributions to 
the TV- fold scattering amplitude matter. These depend 
on the transverse energies and thus require knowledge 
about the transverse motion of the radiating particle in 
the medium. Migdal [|^ was the first to develop a dynam- 
ical description to this aim, employing a two-dimensional 
Fokker-Planck transport equation [ |^, Q for the hard par- 
ton. In the limit of an infinite medium, his well-known 
result shows a characteristic c>c ^/uj low frequency sup- 
pression of dNj/dlogLu compared to the constant depen- 
dence in the Bethe-Heitler limit, i.e., the coherent factor- 
ization limit for a; — > vanishes. This is, however, only a 
very special feature of Migdal's limiting case, where the 
formation length goes to infinity for uj —f 0, but never 
exceeds the (infinite) extension L of the medium. For a 
medium of finite size, the formation length does exceed 
the system size below some critical frequency Wcr, and 
the Bethe-Heitler limit is finite. Migdals -ycJ-dependence 
is hence not valid for cu < Wcr- Moreover, additional 
effects become important for the radiation spectrum at 
lower frequencies [ H . Most notably, this is the transition 
radiation and the Ter-Mikaelian effect [ H] of dielectric 
suppression. 

The renewed interest in the LPM-effect has at least 
two reasons: On the one hand, fourty years after discov- 
ering the theoretical principles [ |lj, pi , the first precision 
measurements of the LPM-effect [0, ^ (and the Ter- 
Mikaelian effect [ §, 1^ ) were made recently by the SLAC- 
146 collaboration. On the other hand, with the advent 
of a new generation of relativistic heavy ion colliders at 
RHIC and LHC, the understanding of the non-abelian 
analogue becomes important. 

In the QED case, the experiment explored relatively 
thin targets with L/lf on the order 0.1 to 10, in which the 
transition between Bethe-Heitler and factorization limit 
occurs. For a quantitative understanding, a realistic the- 
ory has to account for the finite extension of the target, 
the multiple elastic scatterings in the target, multiple 
photon emission, and possibly additional complications 
like radiation off structured targets. 

There are at least three modern approaches, which can 



account in principle for the LPM-efFect in realistic tar- 
gets. They implement the eikonal approximation for the 
radiating hard particle in different ways: i) Blankenbecler 
and Drell [ |l^, [l|] (see also Baier and Katkov [ |l|, |l|, [l4| ) 
started from the solution of the Klein-Gordon equation 
for the charged particle in the presence of an external 
field. The solution was approximated to order 1/ E in a 
high energy expansion in which the eikonal path of the 
radiating particle is recovered, ii) Zakharov [ na pro- 
posed a light-cone path integral formulation which de- 
scribes in a coordinate-space representation the trans- 
verse momentum kicks on the eikonal path, mimicking 
the elastic scatterings by an effective dipole cross section. 
Zakharov's work [ [l^, |l^ provides the most accurate de- 
scriptio n of the measured data so far. As we shall see in 
section [IC , his starting point is closely related to a high 



energy expansion of the Dirac equation in the presence 
of an external potential, iii) A third approach is due to 
R. Baier et al. (BDMPS) [ |l|l, who started from the 
radiation amplitudes for iV-fold scattering in the eikonal 
approximation. For QED, the consistency of their ap- 
proach and the work of Blankenbecler and Drell can be 
checked diagrammatically [ |9| . A discussion of the LPM 
effect also exists for non-equilibrium conditions [ ^ . 

In the QCD case, high pt < 10 GeV jets will be one 
of the new probes of the dense matter produced in rel- 
ativistic heavy ion collisions. Deviations from factorized 
perturbative QCD due to final state rescattering are sen- 
sitive to the density and coupling in the plasma and will 
affect key observables like the high momentum tails of 
single particle spectra. The corresponding study of QCD 
radiative energy loss due to final state rescattering was 
initiated by Gyulassy and Wang [ ^, ^. Recently, it 
has been extended most notably in a series of papers by 
BDMPS [11,11,11,11 using equal time perturbation 
theory. Also, Zakharov has pointed out that his formal- 
ism can be adopted to QCD bremsstrahlung with a dipole 
prescription [ |l^, ^7] , and the equivalence of Zakharov's 
formalism with the work of BDMPS was sketched [ pq . 

All these calculations of the LPM-effect in QCD, how- 
ever, are (1) limited to the cases of infinitely many or 
very few {N < 3) rescatterings of the parton, and (2) 
they mainly focus on the angular integrated energy loss 
dE/dx dL in which the transverse momentum depen- 
dence of the radiation pattern is averaged out. This en- 
ergy loss, however, is not a good observable for QCD, 
because the parton shower is not directly observable due 
to hadronization. In addition, the QCD bremsstrahlung 
of hard jets must compete with the hard radiation associ- 
ated with the jet production. To detect modifications of 
this hard vacuum bremsstrahlung spectrum due to final 
state rescattering requires knowledge about the angular 
distribution of the spectrum. In QCD, the transverse 
momentum dependence of the LPM-effect is thus indis- 
pensable for a quantitative understanding of radiative 
energy loss. 



Recently, Kopeliovich, Schafer and Tarasov (KST) [|27 
have used the Furry approximation of the Dirac equation 
in order to account for the transverse momentum depen- 
dence of the radiation spectrum emitted in a multiple 
scattering process. For the case of QCD bremsstrahlung 
radiation, they translate their QED results into QCD via 
an a posteriori dipole prescription (pTl) . This strategy is 
frequently used not only for the calculation of the ra- 
diative energy loss of high pt partons [ fsl 17, [2^, but 



also for the description of nuclear shadowing in the target 
rest frame [ pSl |2^ and for related problems of diffractive 
dissociation of virtual photons [ ^ . 

As we show in what follows, an extension of the KST- 
formalism has great potential for the calculation of radia- 
tive energy loss in realistic scenarios since (1) it provides a 
smooth interpolation between the cases of infinitely many 
and very few rescatterings and (2) it allows to compare 
QED-inspired calculations of QCD radiative energy loss 
with perturbative QCD results. The present work fo- 
cuses on the general formalism, its region of validity, and 
a qualitative discussion of its generic features. It gives ex- 
pressions which allow for the numerical calculation of the 
radiation spectrum as a function of the medium density 
and extension, but it leaves phenomenological applica- 
tions to further publications. 

Our work is organized as follows. In section O, we 
derive the starting point of our discussion, the radiation 
spectrum (2.22). We discuss how a technical complica- 
tion, the regularization of this spectrum, can be dealt 
with analytically, and we calculate the corresponding in- 



tegrated energy loss. Section [II focuses on limiting cases 
of the general radiation formula (2.22). We derive the 
radiation spectra for A^ = 1, 2 and 3 scatterings and 
we show that these reproduce the Bethe-Heitler and the 
factorization limit. In section [V , the dipole approxima- 
tion of the radiation spectrum ( 2.22 ) is discussed. This 
gives access to true in medium properties of the radiation 
spectrum which cannot be obtained from an expansion to 
finite order in the coupling constant. Finally, we discuss 
in section ^ how this method can be extended to QCD 
and how it compares to results from perturbative QCD. 
Our main results are summarized in the Conclusion. 



II. THE LPM-EFFECT IN QED 

An expansion of the LPM-radiation cross section in or- 
ders of the coupling constant is essentially an expansion 
in the number of elastic scatterings, since each elastic 
scattering Mott cross section is proportional to a^^. In 
contrast, a high energy expansion of the solution of the 
Dirac equation in the presence of an external field takes 
to leading order in 1/E an arbitrary number of elastic 
scatterings into account. In this section, we derive the 
corresponding high energy limit of the QED radiation 
cross section for a hard electron, traversing a medium of 



longitudinal density n(^). The physics contained in our 
main result ( ^.22 ) will be discussed in the following sec- 
tions m - 0. Our approach in [I A parallels to a large 
extent that of Kopeliovich, Schafer and Tarasov (KST) in 
Ref . [ ^ . We present the derivation in full detail to intro- 
duce our notation and to discuss all the approximations 
involved in the calculation. In contrast to Ref. [ E7| , our 



radiation cross section (2.22) contains an e-regularization 
whose treatment is discussed in subsection p^IB| . Then we 
show that the KST-formalism results in an angular inte- 
grated energy loss formula from which modifications to 
Zakharov's formalism [ ll^, O can be obtained. 



A. The KST-formalism: Differential cross section in 
the Furry approximation 

We consider a relativistic electron undergoing multiple 
small-angle scattering in a spatially extended medium, 
described by an external potential [/(x). The angular 
dependence of the radiation spectrum for emitted pho- 
tons carrying away a fraction x of the incident electron 
energy Ei is given by the differential cross section 



d^a 



d{lnx) dpi_ dki_ (27r)4 
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if^\ 



(2.1) 



where k^ and p^ denote the transverse momenta of the 
photon and outgoing electron, respectively. The radia- 
tion amplitude for transversely polarized photons is given 
in terms of the ingoing and outgoing electron wavefunc- 
tions *" and *+, 



Mfi ^ / d*a;^ \x,P2)q.- ee~ 



■c 1^1 ik-x ,j,+ 



*+(x,pi) 



(2.2) 



Here, a. = 7o7, and e^^'^' is the adiabatic switch- 
ing off of the interaction term at large distances. This 
term plays an important role in what follows since the 
e ^ limit does not commute with the longitudinal z- 
integration. The wavefunctions ^* solve the Dirac equa- 
tion in the external potential U{pc): 



d_ 
dt 



C/(x) — m 7o + i a • V 



^{x,pi,2)^0- (2.3) 



we rewrite the Dirac 



Following Kopeliovich et al. [ 
equation 

^a2^-"^' + ^'-2zC/(x)- *(x,pi,2) 

= [-ta- (VC/(x)) - t/2(x)] ^{x,pi,2) ■ (2.4) 

This allows for an expansion to order 1/ E, treating the 
right hand side as a small perturbation, solving the cor- 
responding homogeneous problem and including in a first 



iteration the derivative term VC/(x). The result is the 
Furry approximation [ G^, 0] : 



*+(x,pi) = e-'^i*+^Pi^DiF+(x,pi) 



*f(2;,P2) 



_ —iE2t+ip2Z 



D2F-{^,P2) 



m(Pi) 

V2e; 

V^e; 



(2.5) 
(2.6) 



which is exact to order 0{U/E) and Oil/E'^). The build- 
ing blocks of this solution are the differential operators 
Di and the functions F^ whose properties we explain 
now: 

To determine the functions F^^ we neglect the second 
derivative d"^ jdz^ in (2.4). This is allowed since the lon- 
gitudinal distances in a multiple scattering problem are 
much larger than the transverse ones. The left hand side 



of (2.4) reduces then to a two-dimensional Schrodinger 
equation, with corresponding retarded Green's function 



»;i— + T C/ Z2,r2 

az2 2p 



G{z2,r2;zi,ri\p) 
iS{z2-zi)6{r2-ri), (2.7) 



satisfying G(z2 — ^i, r2; ^^i, ri|p) = <^(r2 — ri) and van- 
ishing for zi > Z2- The functions F^ in ( |2.5| ), ( |2.6| ) are 
solutions of this two-dimensional Schrodinger equation, 

F+(x,pi) = / d2riG'(z,r;z_,ri|pi)F+(ri,z_,pi) , 

^"*(x,P2) = / d'^r2F-* {r2, z+,p2) G{z+,r2; z,r\p2) . 

For very early and very late times, i.e. for far forward 
and far backward longitudinal distances z+ and z- re- 
spectively, they satisfy the boundary conditions 



F+(ri,2:_,pi) =exp<^ ipi±ri 
F"(r2,2:+,p2) =exp<^ ip2±r2 



(2.8) 

(2.9) 

This ensures that the wavefunctions 'i'p approximate 
plane waves at asymptotic distances. The differential op- 
erators Di are obtained by including the derivative term 



2 PI 


2 
2p2 



V U{x) of (2.4) in a first iteration of this solution, 
. a • V a- {pi- npi) 



D, 



k 

'k 



2E^ 
z = n ■ X 



2E, 

; Vi 



IP. 



(2.10) 



We now explain how to obtain from the Furry approxi- 
mation an explicit expression of the differential radiation 



cross section (2.1). We work in the ultrarelativistic limit. 
El ~ pi, E2 ~ P2 = [^ — x)pi, bj = xpi. The lon- 
gitudinal axis is redefined to be parallel to the emitted 
photon k, see Figure El The transverse momenta of the 



ingoing {pi±) and outgoing {p2±) electron for small x are 
therefore 

— 1 1 — X 
Pi± = — k±, P2±=P± k_L, (2.11) 

X X 



q± = P2_L - pi± =p± + kj_ 



(2.12) 




nevJ 



f ^ r;^, = [4 - 4a; + 2x'^] ^-777+ 2mlx^ . (2.17) 
4) We split up the Green's functions 
G(zi,ri;z2,r2) = / dr'G{zi,ri; z' ,r')G{z' ,r'; Z2,r2) 

at longitudinal distances z' in such a way that contribu- 
tions from the amplitude ( p.2|) an d its complex conjugate 
part can be paired, see eg. ( 2.18 ) and Fig. ^. 

After these four steps, the radiation probability can 
be expressed as a Fourier transform over pairs of Green's 
functions on whose transverse coordinates the interaction 



FIG. 1. Choice of the coordinate system in which the new 
longitudinal z-axis is ta ke n alo ng the radiated photon mo- 
mentum, see Eqs. (2.11), (2.12). 



vertex (2.17) is acting: 



For the calculation of Afg, we use the following four 
steps: 



1) We do the time-integral in (2.2) which ensures energy 
conservation, Ei = E2 + lo- 

2) We determine the z-dependent phase factor qz for 



(2_^) by observing that the z-dependent phases in (2.5), 
(2.6) come with the modulus of the spatial momentum 
Pi = \Pi\- Using energy- momentum conservation, one 
finds to leading order 1/E 



q ^ Pi - P2 - k = 



1 



2(l-a;)^i If 



(2.13) 



(|M/ip) = 2Re / dri dr dpdv2 dv[ dr' dp' dr'^ / dz 

00 
Z 

^ T^n y2{G{z+,V2]z', p\p2) G* {z+,y'^;z', y'\p2)) 

ipi (1 — a;)^ 

xf_, f;, (G(z', p; z, r\p2) G*(z', r'; z, p'\p,)) 

X (G(z, r; z_, ri|pi) G*(z, p'; z^,v[\p,)) . (2.18) 

To turn this into an explicit expression, it is necessary to 
i) calculate the averages (...) over the scattering centers 
in the medium and ii) evaluate the Green's functions G 
which are given by path integrals 



The inverse momentum transfer 1/g is often refered to as 
photon formation length If. However, contributions from 



Giz',r';z,r\p) = 



transverse energies k^/2w dominate the z-dependent / 'Dr{^) exp < / d^ 



phase factor oi Ma, as will be seen below. Hence, If gives 
only a rough upper estimate for the longitudinal size over 
which interferenc e eff ects suppress photon emission. 
3) We absorb in (2_^) the spinor structure of the wave- 



*P.2 



(0-*C/(r(0,e) 



(2.19) 



functions ^J, ^^ in the interaction vertex 



Tr = Vl — xu*{p2) D2 a ■ eDi u(pi) . 



(2.14) 



Depending on the electron spin of the ingoing (Ag — zt^) 
and outgoing (Ae' = i^) electron and the photon helicity 
(A^ = ±1), this vertex takes the form 

Tr (Ae = Ae' , A^) = —i A^ [A-y (2 — x) -I- 2 Ae x] 



with boundary conditions r(z) = r, r(z') = r'. We show 
in Appendix H that the in medium averages involved in 
( ^.18 ) result in a factor which depends on the cross sec- 
tion a (p(C)) times the density n(^) of scattering centers 
as follows 



e^plz d^ [U{r{0,0-Uir'iO,0] 



= exp 



rfes(e,[r(c)-r'(c)]; 



d d 

X I t; * A-y — — 

Fr (Ae = - Ae' , A^) = 2 me a; A.^ Sx^ ^K 



(2.15) 
(2.16) 



T.{i,p{0) = i^n{(,)a{p{0) 



(2.20) 



(2.21) 



where the spin-flip contribution ( p. 16 ) is only non- 
vanishing for A-y = 2 Ae . The differential operators in 
( p. 14 ) act to the right on the two different transverse 
components of r = {rx,ry). In what follows, we are only 
interested in the spin- and helicity-averaged combination 
Tr F*, which takes the simple form 



An explicit expression for cr(p(^)), supporting its inter- 
pretation as a Mott cross section, is also derived in Ap- 
pendix Rl Most importantly, the resulting expression 
^(Ci PiiT) depends only on the relative distance p(^) be- 
tween the paths r(^) and r'(^). As we explain in Ap- 
pendix H, this makes it possible to carry out a large num- 



ber of the integrals and path- integrals in (2. IS) analyti- 
cally. One obtains (see Appendix M for further details) 



d^c 



(i(lnx) dY>± dkj 



/' 



= Cprc Re / dp2 dpi 



I dz I dz'cxp{— ig(z' — z) — e(|z| + |z'|)} 

Z— z 

(Z 2 + 

- j di^{^,Xpi)- j d^n^,xp2) 



1-X 



X exp <—ix\ p^ 

d d 
dpi dp2 



fc 



• P2 - ik_L • Pi 

/C(z',P2;^,Pi|a*)- (2.22) 



This differs from the resuh in [ E^ by including the reg- 
ularization e"'^'-'^'"'''^ '^ The physical radiation spectrum 
is obtained by removing e -^ after doing the z and 
^'-integrations and taking z_ -^ — oo and z+ — > oo. To 
simplify the notation, we define 

li = Ei{\-x)x, (2.23) 

which plays the role of mass, and the factor 

2x2 



P" = (2^ El (1 - x)2 



(2.24) 



which accounts for a combination of recurring kinemat- 
ical prefactors. The coefficients of the spin-flip (^gf) 
and n on-flip (ftif ) contributions of the interaction vertex 
(^.17|) read 



5nt 



4 - 4a: -I- 2a:^ 
4a;2 



5sf 



(2.25) 



The remaining path integral K, in (2.22) is given by 
K.{z',Vc{z')\z,Yc{z)\ti) 



= / Vyc exp 



d£. 



-Yl+iY.{S„xVc 



(2.26) 



In this path integral, S(^, a; Tc) plays the role of an imag- 
inary potential, while its physical interpretation is up to 
a minus sign that of an clastic Mott cross se ction times 
the density of scattering centers, cf. equation (A12). De- 
pending on the context, we shall hence refer to S as a 
potential or a cross section. 

The following sections are devoted to a study of the 
differential cross section ( ^.22| ). We focus on analytically 
accessible limiting cases which illustrate the physics con- 
tained in (2.22) and we discuss approximations which 
allow for its numerical analysis. 



B. Removing the regularization e 







The regularization prescription in (2.22) cannot be ne- 
glected since the e ^ limit does not commute with 



the z- and z'-integrations. This complicates practical 



applications of (2.22): even if the Green's function K, is 
known explicitly in some approximation (see sections III 



and IV below), the expression (2.22) is not suited for nu- 
merical calculations since one cannot control numerically 
the e — > 0- limit after carrying out the z-integrations. 

To solve this problem, we determine the e — > limit of 
the radiation spectrum ( ^.22| ) analytically. We consider 
a medium of arbitrary but finite longitudinal extension 
which is positioned along the longitudinal axis between 
and the finite distance L. The z- and z'-integrations of 
(p. 22) can be split into six parts. 



Z+ 2+ 



L 



2+ 



z z- z- L 

L L L Z-y ^+ ^ + 



(2.27) 



L 



L z 



To shorten the following calculations, we omit the spin- 
fiip contribution g^i to the radiation spectrum. This con- 
tribution is negligible in the relativistic limit rrie <C Ei 
on which we focus in the following. [If needed, the g^i- 
contributions to all our results can be recovered without 
additional technical difficulties.] 

The differential cross section receives six contributions 



which we label according to (2.27) in an obvious way 



dV 



— L/pio 5nf 



d{\nx) dp^ dkj 

X (/1+/2 + /3 + /4 + /5+/6) 



(2.28) 



Since the medium does not extend outside the inter- 
val [0, L], we can replace the in- medium propagator 
]C(z', P2] z, p il/i) outside [0, L] by the free propagator /Cqj 
given in (3.3) below. E.g., for Z2 > L > zi, we can write 



IC{z2,P2;zi,pi\n) = / drJCo{z2,P2;L,r\fi) 

xlC{L,r;zi,pi\fi). (2.29) 

This allows to do all "infinite" integrals J^ and J^^ in 
( p.27D analytically. One is left with the "finite" inte- 
grals /„ for which the e ^ limit commutes with the 
z-integration. Doing the z-integrations first, then tak- 
ing z_ — !■ — CX3 and z+ — » +00 respectively, and finally 
adiabatically switching off the regularization, e ^ 0, we 
find 



..,..^)/.„-r— ^^-^^ 



(2.30) 



I2 = -Re / dz'e"'«^ / dre 




- I S(C.a:r)-Jxp2i-r ^ 



X /dfe"Pi^-^a;pi^-— /C(z',r;0,f|Ai), (2.31) 

h = (a;pi±) • (a;p2±) Re I dn dra e-^^P^^""^ 

xe-Pi^"-i/c(L,r2;0,ri|/i)Z3e-'^«, (2.32) 



h^Re f dz f dz' fdridr2e-"P 






2i-r2 I2:pii-ri 



(2.33) 



1j 



X J dre-'^P'^-' xp2± ■ ^ICiL,r;z,r\fi) , (2.34) 



/e = (x^ p2^ 



dre Jo Z, 



6 • 



(2.35) 



Here, the factors Zi contain the explicit solutions of those 
z-integrals which extend to z+ or z- respectively. Using 
the notational shorthands 



Qi = 

we obtain 



9 9 



n2 /'r,2 



q = 



2fi 

x^ (pL + ml 
2iJi 



(p i± + o 

2fi 



(2.36) 

(2.37) 



u u 



1 



2Q2' 

and analogously, 

1 
Zo^ — , 
Qi' 

1 

' Q2' 



23 

26 



1 



Q1Q2 

1 



2Q2- 



(2.38) 

(2.39) 
(2.40) 



The differential radiation spectrum (2.28) given by the 
equations (2.3C) - (2.40) is suitable for a numerical eval- 
uation since the e — > 0-limit is taken. Previous numer- 
ical investigations [ M of ( p. 22 ) were based on a non- 
regularized expression which did not include the term 
g-e(|z|+|z I) jj^ ^Yie integrand. The difference to the above 
results is found by determining the above z-integrals 
without regularization: 



u u 

Zl''' = Re f dz f dz'e~"^^^''-''> 



= 2 2i (l-cos{Qiz_}) , (2.41) 

Z^"^ = Z2 (1 - exp {z Qi z_ }) , (2.42) 
Zr = 23 (l-exp{iQiz_}) 

X (1 - exp {i Q2 (2;+ - i)}) , (2.43) 

Zr = 25 (1 - exp{z Q2 {z+ - L)}) , (2.44) 

Z^''- = 2Z6 (l-cos{g2(z+-L)}) . (2.45) 

The differential cross-section ( [j.22 ) can be related to the 
unregularized one cr("'") : 

d{lnx) dp± d]i± d{lnx) dp± dk± 



C^ 



dpcxpl- I](^, xp)-ixq±- p 



{x^pIj 



{x-pl^) 



Qi 



6° 



[z_,z+]. (2.46) 



Here, the correction term 0°q^^[2;_, z+] summarizes those 
oscillating contributions which vanish if one averages 
suitably in a numerical calculation over the boundaries 
Z- and z+ of the longitudinal integration. In addition, 
there are two terms which do not depend on the bound- 
aries z_, z+. They stem from the terms ZJ"" and Z^'' 
which differ by an overall factor 2 from the regularized 
expressions Zi , Zg . These terms provide the proper sub- 
traction of the Weizsacker- Williams fields of the ingoing 
and outgoing electron. As will become clear below, they 
are n on-negligible, and ensure e.g. the consistency of 
(^.22 ) with the Bethe-Heitler spectrum in the appropri- 
ate limiting cases. 



C. Integrated energy loss cross section 



The total energy loss cross section is calculated by in- 



tegrating (2.22) over transverse momenta. To this aim. 



we rewrite the differential cross section as a function of 
the relative momentum transfer q_L = k_L -I- pj_ and the 
transverse momentum k^ of the emitted photon. Since 
there is no kinematical boundary on q^, the integration 
over the transverse momentum of the electron is simple: 



dh 



d(lna:) dkj 



(2^) 



C^ 



Re / dpi dz dz'e-'«(^'-^) 




dpi dp2 



Z 

/ dCS](C,xpi) -ik^ ■ pi 



IC{z',p2 = 0;z,pi\fi). (2.47) 



The integration over the photon transverse momentum 
is more complicated due to the finite kinematical bound- 
aries for the photon, |kj^| < lu — x Ei: 



da 



d{lnx) 



rfk 



dV 







d(hix) dkj 



a. 



cm 



2 Re 



(27r)2 El (1 - a;)2 



dp_ 



2 TT U) Jl{ul pi) 



Pi 






xexpl- S(C,a;pi) 
xK{z',p2 = 0;z,pi|^) 



9 a 

dpi dp2 



(2.48) 



It is interesting to compare this energy loss formula to 
that derived by Zakharov [ 15, 17| : 



^^(zak) 

d(lna;) 



2 Re 



El (1 - xY 






9pi dp-2 
X [/C(z',0;z,0|/^) -/Co(z',0;z,0|^)] . (2.49) 



dpi dp2 



Here, /Co is the free, non-interacting path integral, given 
explicitly in (3.3) below. If we had ignored the kine- 
matical boundary, extending the fcj^-integration in (2.4S) 
up to infinity, we would have regained the expression of 

Zakharov except for the term proportional to /Cq. 

We note that Zakharov's arguments leadi ng to ( 2.49| ) 
are very different from our derivation of (2.48). The 
transverse momentum dependence of the radiation enters 
in no intermediate step of his calculation [ |lj. Also, the 
/Co-term in ( 2.49| ) does not result from a derivation: it is 
rather subtracted a posteriori as a 'renormalization pre- 
scription' in order to cancel a singularity in the integrand 
of ( |2.49| ) for small {z' — z). In t his se nse, our derivation 
differs from Zakharov's result ( 2.49| ) by: i) taking the 
finite phase space of the emitted photon properly into 
account, ii) containing a proper e-regularization of the 
radiation amplitude and iii) not employing a subtraction 
of a singular contribution a posteriori. 



III. THE LOW OPACITY EXPANSION FOR 
THIN TARGETS 

Information about the target medium enters the radi- 
ation cross section ( 2.22] ) via the product E(^,p), which 
measures the elastic cross section times the density of 
scattering centers in th e me dium. Since S(^,p) oc a^.^^ 
n(^), an expansion of ( 2. 22 ) in powers of a^^ is an ex- 
pansion in powers of the opacity Tucff', where Ccff is the 
effective elastic single scattering cross section and 



r 



d^niO. 



(3.1) 



T CTeff measures the average number of scatterings for an 
electron traversing a med ium of length L. Since the elas- 
tic cross section in ( p. 22 ) is not a geometrical quantity 
but depends on the integration variables, our expansion 
of ( 2.2^ ) will be formally in powers of a^^T, and we 
shall refer by a slight abuse of language to T as opac- 
ity. Of course, after all integrals are done, each power 
of al^ T will be accompanied by a power of the e lastic 
cross section. The A^-th term in this expansion of ( 2.22| ) 



is of order aom T^ and corresponds to the multiple 
scattering off exactly N external potentials. An expan- 
sion in the opacity T is thus an expansion in the number 
of rescatterings. 

Below, we study this low opacity expansion of the radi- 
ation spectrum up to third order. This leads for equation 
( [2.22 ) to a number of consistency checks which any radia- 
tion spectrum including in medium effects should satisfy. 
Moreover, this provides the basis for a discussion of the 
corresponding QCD radiation spectrum in section M. 



A. Bethe-Heitler cross section as a low density limit 



The derivation of the radiation spectrum ( p.22[ ) relies 
on the approximation that the distribution of scattering 
centers in the medium can be described by the average 
( ^.20 ) . This uses explicitly that the size of the medium is 
much larger than the Debye radius of a single scattering 
potential. A priori, it may hence seem unclear to what 
extent one can still recover from ( ^.22 ) the correct radi- 
ation spectrum for a single scattering process where the 
extension of the potential is the extension of the target. 
However, in the low opacity limit, when the distance be- 
tween scattering centers is much larger than the photon 
formation length, multiple scattering should not affect 
the radiation pattern. The radiation spectrum should 
converge to the single scattering Bethe-Heitler cross sec- 
tion times the opacity factor T. Deriving this is an im- 
portant consistency check for our formalism. 

Our expansion of the integrand of the cross section 
( ^.22 ) in powers of S(C,p) oc a^.^ "(C): ^ses the corre- 



sponding expansion of the path integral 
/C(z', p2;z, pi) = K.q{z\ P2; z, pi) 

I 

Z 

- dz dplCo{z\p2;z,p)'E{z,xp) ICo{z,p;z,pi) 

Z 

z' z' 

+ dzi / dz2 / dfidf2JCo{z\p2;z2,f2)T.{z2jXf2) 

Z Zi 

xJC{z2,f2;zi,ri)j:{zi,xfi)JCo{zi,ri;z,pi), (3.2) 



where from now on, we suppress the expHcit fi- 
dependence in /C and in the corresponding free Green's 
function ICq, 



ICq{z',P2;z,pi) 



M 



27ri(z'-z) 



exp • 



ifj. (pi - Pa) 

2(z'-z) 



(3.3) 



From the two potentials in the third hne of ( 2.22) and 
from the propagator /C, the radiation spectrum ( 2.22| ) 
receives three contributions to order 0{al^^T): 



d^a 



dilnx) dpj^ dkj 



0{alr^ T) 



Cp,,gni ih+h+h) ; (3.4) 



2+ 2+ 



z_ z 

z 

xe"P-P2-..Pi^PiyS(^,:cpi)Ki(pi,p2), (3.5) 

z_ z 

xe^-P-P^-"Pi^PiyE(^,xpi)Ki(pi,p2), (3.6) 

z' 

Z+ 2 + 

h = -Rc /" dz /" dz' e-"?(^'-^) /" rfp2 dpi 



xe 



lXp2J_-p2-2 3:pii-pi 



'«2(Pl,P2)- 



(3.7) 



Here, ki and K2 denote the derivatives of th e ze roth and 
first order contributions to the propagator ( p.2| ), 



Ki(Pi,P2) = -R — ■ ^ — ^o(z',P2;z, Pi) • 
9pi 9p2 



(3.8) 



K2(Pi,P2) == -Q -Q — / iz \ dpK.Q{z',p2]z,p) 

xj:[z,xp)K.o{z,p;z,pi). (3.9) 

We consider a medium of finite longitudinal extension L 
and homogeneous density n{^) = no- Then, all integrals 



in (3.5) - ([3.7|) can be done analytically. Taking the limit 



e -^- after the z-integrations, we arrive at 

h - -^^0^^ ^l.(qi), 

^ ^ x2p?^(27r)2 



I3 = -L no 



2Ql a;2 
Qi Q2 x^ 



S(q±). 



(3.10) 
(3.11) 
(3.12) 



The result is particularly transparent in the relativistic 
limit, when all mass dependencies can be neglected and 



g « 0. For a medium characterized by a set of Yukawa 
potentials with electric charge Z e and Debye screening 
mass M (see Appendix ^ for a derivation), 



- S(qx) 



4Z2aL 



q-L#o 



[A-P+cil) 



2 ^2 



(3.13) 



den otes the elastic scatterin g cro ss secti on with qx given 
by (|2.12| ). The three terms (|3.10| )-( |3.12| ) then combine to 



d^c 



d(lna::) dp^ rfkj 



0(q3„ t) 

16a;2 



qi 



(2 7r)2 (M2 + q^)2 ki(ki-.Tqx)- 



■ r. 



(3.14) 



This is the high-energy limit of the Bethe-Heitler radi- 
ation cross section for scattering off a Yukawa potential 
times the opacity factor T = Luq. The exact result 
[see e.g. equation (5.150) of Ref. [ |3^ in the coordi- 
nate system Fig. p] used here.] differs from the above 
approximation only by terms of order 1 -f 0{q^ /E^) and 
1 4- 0{x^). The main properties of the QED radiation 



spectrum are seen clearly from (3.14). The spectrum van- 
ishes for vanishing momentum transfer q±^ — > 0, it shows 
the x^-dependence characteristic for the QED radiation 
spectrum peaking at forward rapidity, and it shows the 
correct 1/k^ (fej_ — xq±)'^ dependence. 

We mention as an aside that the evaluation of the cross 



section (3.4) without regularization prescription leads to 
terms /i and I2 which are a factor 2 larger while the inter- 
ference term remains unchanged. The e-regularization of 
( ^.22 ) is thus necessary to regain the Bethe-Heitler spec- 
trum. 



B. Low N multiple scatterings 

In complete analogy to the derivation of th e Bet he- 
Heitler limit, one can expand the cross section ( 2.22 ) to 
second order. All contributions to the radiation cross 
section depend now on the product of two elastic Mott 
cross sections S, whose combined momentum transfers 
sum up to qx. For this we introduce the shorthand 



dV^{q^)= / dqiidq2±S(qix)S(q2±) 

X^<5(2)(qi-qii-q2±). (3.15) 



Also, we introduce shorthands for the transverse mo- 
menta and the corresponding transverse energies, 

ui = a;pi± , U2 = xp2i_ , 

u^^x (pii + qii) = X (p2_L - q2±) , 
222 

(O — "1 ^ _ "m ^ _ "2 
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The second order contribution of the radiation spectrum 
( ^.22 ) consists of six terms. After integrating out aU 
transverse coordinates, it takes the form 



d\ 



d{hix) dp± d]s.± 



o(ccl^r^) 



Cprc5nf / dV'^{qA_) 



2^(2) 



2 1 1 



<22 



(2) 



2^(2) 



2^223 



r(2) 



-U„ 



U2Z, 



(2) 



U„ 



UiZ, 



(2) 



Ui • \l2Z\ 

]. (3.16) 



The first three terms stem from the expansion of the ex- 
ponential term exp \— ^^ S — JJ^,^ E I in (2.22) to sec- 
ond order in E, the fourth term is from the second order 
expansion of the Green's function /C, and the remain- 
ing two terms are contributions from the first order in K, 
times the first order of the exponential term. The vari- 

("2) 

ables Zl stand for the remaining longitudinal integrals 
over phase factors. In general, they are involved, e.g. 



2+ 
Zf^ ^Re I dz 



^+ 



d/e-^(l^l+l^'l) dzi 



no 



dz' 



Xe 



-iQ2(z' -Z2)-iQm(z2-zi)-iQi(zi-z) 



n{z2) 
no 

(3.17) 



For a medium of homogeneous density no and finite 
length L, they reduce to very simple forms: 



r(2) 



r(2) _ 



r(2) 



2gf' 

cos(LQ,„)-l + 5i'Qm 



L2 



2Qr 



,(2) _ COs(LQ,n) - 1 



r(2) 



Q1Q2QI 

cos{LQm) 



SL^Ql 



Q2QI 



(2) _ cos(L Q,„) - 1 + iL2 g^ 

^6 — 



QiQl 



(3.18) 
(3.19) 
(3.20) 

(3.21) 
(3.22) 
(3.23) 



From this, we read off simple limiting cases: In the limit 
of a very thin target of fixed opacity, we can move the two 
scattering centers so close together that the photon can- 
not resolve them. The spectrum is then indistinguishable 
from a single scattering process 



d^a 

lim 

L^o d(ln.x) dpj^ dkj^ 

^2 



r fixed 



0(al^T^) 



L-prc 5nf ■ 



Ul 



U2 



4 \Qi Q2 



dV^qj 



(3.24) 



This is the coherent factorization limit that corresponds 
to one single effective Bethe-Heitler radiation spectrum 



which depends only on the initial and final momenta. 
The characterisitc momentum dependence of the Bethe- 



Heitler radiation spectrum (3.14) is now combined with 
the convolution of two Mott cross sections ( |3.15 ). The 
first L-dependent correction to this factorization limit is 
proportional to L^ and takes the form 



Cpregm'^L^ JdV^{q^)Ql 



ul 



Ul Uz 



U2 



u„ 



Ul 



Qm Ql Q2 Qm Q2 Ql Qr. 



(3.25) 



which is mainly of formal interest since the q^-integral 
diverges logarithmically. 

In the opposite limit, i — )■ 00, we can study for fixed 
opacity T = Lhq the case of two well-separated scatter- 
ing centers: 



lim 

L^oo d{lnx) dp± dkj 



T' 



U]_ 



Ur, 



T fixed 

o(Q=,„r2) 
2 



a 



u„ 



'pre5nf / dV^{q±) 

(3.26) 



U2 

Q2 



This is the incoherent Bethe-Heitler limit in which the ra- 
diation spectrum is the sum of two independent Bethe- 
Heitler contributions associated with the scattering off 
the first and the second external potential with momen- 
tum transfers qi^ and q2±, respectively. It provides 
the starting point for a comparison of the corresponding 
QCD radiation spectrum with a full N = 2 pcrturbative 
calculation in section ^. 

The pattern emerging in this expansion to second or- 
der in the opacity is fully confirmed in the case of A^ = 3 
scattering centers. In Appendix y, we derive the cor- 
responding radiation cross section and check the Bethe- 
Heitler and factorization limit. 



IV. DIPOLE APPROXIMATION FOR TARGETS 
OF ARBITRARY EXTENSION 



To calculate the radiation spectrum ( 2.22 ) to all or- 
ders in opacity, an approximation scheme for the path 
integral /C is needed. The low opacity expansion, stud- 
ied in section III, becomes rapidly more complicated for 
increasing number of scattering centers. It is only use- 
ful for the description of ultrathin media where very few 
elastic scatterings have to be taken into account. For 
the generic case of a hard projectile particle undergoing 
many small angle scatterings, the dipole approximation 
of /C is standard [ |l|, 0, |l|, |2|, H, |3|] . 

This dipole approximation is based on the observation 
that for small transverse distances p — \p\, the cross sec- 



tion a{p) of the Yukawa potential (All) has a leading 
quadratic dependence: [ 13 



C{p)^An{Za,^f 



^""^^ITp 



1-27 



(4.1) 
(4.2) 



where 7 — 0.577 denotes Eulers constant. The main con- 
tribution to the radiation spectrum (2.22) comes from 
small values of p, where C{px) shows only a slow log- 
arithmic dependence on p and can be approximated by 
a constant, C = C{pcsx). For a quadratic dependence 
cr{p) = C p^ , where this logarithmic dependence is weak 
enough to be neglected, the path integral K. in ( 2.26D is 
that of a harmonic oscillator f [151 



JCos7.{z2,'r2\Zi,Yi 



fxQ 



X exp ■ 



27ri sin(riAz) 
ipfl [{rj + rj) cos(r2Az) - 2 ri • rs] 
2 sin(f7Az) 



(4.3) 



with the oscillator frequency 



n = 



1 — i nC x^ 



V2 



M 



(4.4) 



In the following, we exploit the consequences of this ap- 
proximation. 



These terms show a characteristic dependence on the cou- 
pling constant: in the dipole approximation, C is con- 
stant and measures the size of the elastic single scatter- 



ing Mott cross section, see Eq. (4.2). Thus, C and VP 
are proportional to aoim ^'^d the three terms I\, Iq and 
/s are non- analytic in the coupling constant. For small 
ttpffj, their leading dependence is a~^ x exp {©(a^Tm)}. 
We show in the next subsection that these contributions 
combine to a radiation cross section with a leading de- 
pendence acm X 



exp {©(a^jjj'^)}. The factor 



exp 



qi 



2nCL 



(4.9) 



responsible for the remaining non-analyticity of the radi- 
ation cross section has an intuitive physical interpretation 
as Moliere q^ -broadening. This broadening is propor- 
tional to the density of the medium no, the path length 
L inside the medium and the probability C that a scatter- 
ing center in the medium interacts with the hard electron. 
As expected from a particle undergoing random motion 
in the q_L-plane, the accumulated average q^^ grows pro- 
portional to L. The corresponding contribution (4.9) is 
seen explicitly in /i and /g, but it can also be recovered 
from I3 in an expansion to lowest order in a^^, when 



i pfl sin(f7 L) = nC Lx^ + 0{a 



^ ) 

cm/ 



(4.10) 



A. The general expression: Non-analyticity in the 
coupling constant 



We now turn to the remaining three contributions of the 
radiation cross section (4.5), which arc given in terms of 
integrals over the longitudinal extension of the medium: 



In section [IB, we have determined the six contribu- 
tions (2.3C) - ( 2.35| ) to the radiation spectrum ( 2.221 ). In 
the dipole approximation, the transverse integrations in 



(2.30) - (2.35) reduce to Gaussian integrals and can be 
done analytically. Here we discuss the resulting expres- 
sion 



dV 



(i(lnx) dpj^ dkj 



— C'pre ffnf 2^ I 



(osc) 



(4.5) 



i=i 



obtained from (2.28). For three of the six terms R 
all integrals can be done analytically: 



(osc) 



r(osc) 



r(osc) 



1 p2^ 2^ 



2 Q{ nCL 

1 p2^ 2^ 



2 Ql nCL 



exp 



exp 



qi 



2nCL 



qi 



2nCL 



(4.6) 
(4.7) 



^(osc)^ X Pl^^p^i. ^^^^,^, 



Q1Q2 

X exp 



ix'^ (1 — cos(r2L)) ^^2 



2pn sin(f7L) 



(pL 



Pi±. 



2i7r 



p^ sin(il L) 



exp 



9 ? 
IX q]^ 



2pn sm{nL) 



(4.8) 



r(osc; 



— T^ aze 



CI 



,(1). 



,(12) 



Pl± • P2± 



exp 



r(osc) 



4 Co 



5^ ^ p^^ - 2 pii • P2_L + &: 



(2) 







,(2) 



exp 



r(osc) 



P2± + Og 

2 



(12) 



Pl± P2J 



— X 

4C; 



4S?, 



-2pi 



P2_L 



L 

Re / dz 



dz'e-"^^^'-'^ 



pL 







,(0) 



,(1) 



Pij 



,(12) 



Pl± ■ P2J 



exp 



4C4 



b['>pl^-2p. 



P2_L 



+ bip 


vli_ 




CI 


X 


+ af 


P2± 


+ bf 


pL 



(4.11) 



(4.12) 



(4.13) 



To make the structure of these terms transparent, we 
have introduced several shorthands a^, bi and Ci. The 
factors bi and Ci show simple limiting behaviour: 
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6, = 1 + OialJ , 
AC, ^2nCLx^ + 0{ai^) . 



(4.14) 
(4.15) 



This allows us to recover the Moliere qj_-broadening term 
(4.9) to leading order in the coupling constant in all six 
contributions n°^'^' of the radiation spectrum (4.5). 
The explicit form of the functions Ci are: 



C2 = ^-^n{n{L- z) cos[rj z] + sin[f7 z]} , (4.16) 

C5 = ^-^n{nz cos[n{L - z)] + sm[n{L - z)]} , (4.17) 

Ci^^-^n {~n^ z{L~- z) sin[r2(z' - z)\ + sin[f7(z' - z)\ 



+r2 (L - z' + z) cos[r2(z' - z)]} , 



(4.18) 



from which the limit (4.15) is recovered easily. The func- 
tions hi are given by 



^(2 



^(2 



^(2 



= cos[r2 z] — fi (L — z) sin[f7 z] , 

= cos[i7 z] , 

= cos[r2 (L — z)] , 

= cos[J7 (L — zy^ — ^z sin[rj (L — z)] , 

= cos[r2 (z' - z)] - f7 (L - z') sin[r2 (z' - z)] , (4.23) 

= cos[rj (z' - z)] - f7 z sin[r2 (z' - z)] . (4.24) 



(4.19) 
(4.20) 
(4.21) 
(4.22) 



from which the limiting form (4.14) is also clear. Finally, 
the prefactors of the integrands in ( 4.11 )-( 4.13 ) are: 



2 2 

2 



r22(L_2), 



(12) _ M 



f7 sin[r2 z] 



,(2) 



= ^f^^z. 



(12) A* p, . r^,^ 

a\ = —\l sin[S2 (L 



x^ 






,2n4 



/^2^ 

4 

/i2r)4 

4 

,2o4 



(L-z')"c + 



f7 



x((L-z')c+l)} 



(4.25) 
(4.26) 
(4.27) 
(4.28) 

(4.29) 
(4.30) 



where we have used c = cos[r2(z' — z)] and s = sin[r2(z' — 
z)]. The expressions given here allow for a numerical cal- 
culation of the radiation spectrum, once the extension 
L of the medium, its density np, and the measure C of 
the elastic cross section of a single scatterer are given. 
Based on these expressions, we discuss in the next sub- 
sections i) an analytically accessible limiting case of the 
spectrum (4^) and ii) a Bethe-Heitler limit which allows 
to determine the constant C phenomenologically. 



B. Moliere limit: Radiation in the Gaussian small 
angle multiple scattering regime 

We now turn to a limiting case of the radiation spec- 
trum (2.22) in which a characteristic in-medium effect 
can be isolated in a simple analytical expression. This 
limit focusses on the kinematical region 



(7_L < fc± < fc < ^ , 

9i> 1, 
|f^L|<l. 



(4.31) 
(4.32) 

(4.33) 
(4.34) 



The first condition ( 4.31 ) is satisfied in the high energy 
limit when x is small and the photon is radiated under a 
small angle with respect to the beam. It is characteristic 
for relativistic kinematics that the transverse momentum 
freed by multiple interactions can be significantly larger 
than the total transverse momentum q^ transfered by 
the medium, i.e., p^ « — k^ and g^ = |p^ -|- k^| ^ k^. 
According to ( 4.32| ) and ( 4.33| ), the length of the medium 
has to exceed one formation length significantly, while the 
last condition |riL| <C 1 can be realized e.g. by choosing 
a target of sufficiently low density. We note that (4.33) 
is by far the most stringent condition. 

As a consequences of ( [1.32 ) and (4.3E), the terms 
cos{QL), sm{Q L) are rapidly oscillating as function of 
L. If we assume a target with varying extension L, e.g. 
due to an unpolished surface, then contributions to the 



radiation spectrum (4.5) proportional to these oscillat- 
ing terms will be averaged out by the experiment. An 



expansion of the terms in (4.5) in powers of the coupling 
constant leads then to 
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e 2 ,. c t 



= 
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nCL Q2 
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27r 



27r 47r 1 

Q2 + Q?I^ 

47r 1 



2 Q2 ^(jL 



Q2 L2 Q2 



cxp 



e 2„ci ^ 



2nCL 



(4.35) 

(4.36) 

(4.37) 
(4.38) 

(4.39) 

(4.40) 



Here, we have used (4.32) to neglect the qj^-dependent 
terms except for the leading 0(1/q;2^j) Mohere factor 
(4.9), and we have approximated pi± « p2_L ~ P*_l- 

2 

Corrections to these terms are of order 0(|^) and 
0(1/(2.2 Q2)2)^ Neglecting 0{1/{L^ Q'^))-contTihutioiis, 
the radiation spectrum (4.5) takes the simple form 
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d^cr 



d{\nx) dp±^ dkj 



6 

— r* ^ Y^ 7-(osc) 

— L^prc i/nf / ^ -(j 



ao 



exp 



2nCL 



(4.41) 



This Moliere limit is the main resuh of the present sub- 
section. In accordance with the Bethe-Heitler single scat- 
tering cross section (5.14), we still find the characteristic 



X rapidity dependence, as well as the l/k]_ fall off of 
the spectrum in the kinematical regime q± <C kj_. In 
contrast to the single scattering Bethe-Heitler cross sec- 
tion, however, the spectrum now peaks at vanishing to- 
tal momentum transfer q± = rather than to vanish for 
qx — * 0. The reason is that the radiation from many 
small angle scatterings adds up to a finite contribution, 
while the sum ^,- q^j^ of many small random momentum 
transfers q^^ peaks at zero. 



Also the dependence of ( 4.41 ) on the coupling constant 
is easy to understand. The cross section is proportional 
to aom, since one photon is radiated off. The broadening 
of the q_L- distribution is determined by the probability 
that the hard electron undergoes an interaction. This 
probability is given by the elastic single scattering Mott 
cross section, and hence the exponent shows a character- 
istic a~f^-dependence. 



The radiation cross section (4.5) obtained in the dipole 
approximation, is clearly more complex than the limiting 
case (4.41). The latter, however, illustrates most clearly 
that generic in medium effects are retained in the dipole 
approximation which cannot be obtained from a calcula- 
tion to fixed order in the coupling constant. 



C. Fixing the dipole parameter C in the 

Bethe-Heitler limit: validity of the dipole 

approximation 

In the dipole approximation, the expansion of the radi- 
ation cross section ( 2.22 ) in powers of the coupling con- 
stant does not converge. The spectra for N = 1,2,3 
scatterings derived in section [II, are proportional to the 
N-th power of the elastic scattering cross section S(q±) 
which is essentially the Fourier transform of cr(p). While 



the Fourier transform of the analytic expression (All) for 
(j{p) is well-defined, the Fourier transform of its dipole 
approximation cr(p) — C p^ diverges. 

This failure of the dipole approximation to allow for an 



expansion of (2.22) in powers of the opacity T does not 
affect its validity for the calculation of medium effects as 
given e.g. in the last subsections. This can be seen e.g. 
from the factor exp {— J E(^, a; p)} in ( p. 22 ). Its Fourier 
transform is well-defined and it is well approximated in 
the dipole approximation (7{p) = C p^ by a Gaussian of 
appropriate width. It is only the expansion of this factor 
in Qfcm whose Fourier transform does not receive the main 



contributions from small values around p — Q. This is the 
reason why the dipole approximation of <t{p) for small 
values of p cannot be combined without difficulties with 
an expansion in the coupling constant. 

For the q^-integrated spectrum (2.47), this technical 
difficulty docs not exist. The reason is that in the in- 



tegrand of (2.47), the first order term E(^,xpi) comes 
multiplied by K.{z' , 0; z, pi|/J.) which is to lowest order in 
the coupling constant a Gaussian in pi . This ensures that 



the main contribution to the integral (2.47) comes from 
small values of p where the dipole approximation is valid. 
The corresponding Bethe-Heitler spectrum in which the 
dipole approximation is employed after integrating out 
the qj^-dependence, reads: 



dh 



d{\j\x) c?ki 



dipole 



Oia^^^^T) 



'-*cni 2 



C {4-4a; + 2a;2} 
k^ -(- (r?T,e xY 



{\^l + {m,xyy 



T. (4.42) 



In the calculation of ( 4.42| ) we have no t neg lected the 
g-dependence of the radiation spectrum ( 2.22|) . The con- 
sequence is the apprearance of the very small term rUe x 
which regulates the -^ singularity. This allows to calcu- 
late the kj^-integrated Bethe-Heitler energy loss formula 
in the dipole approximation 



da 



d{\nx) 



dipole 



0(a?,„ n) 



1 2 2 

1 — m„x 



,C 



2)11 ml 



{4-4.T + 2x2} 



iuj^ + 3w^ rngX^ -I- TOgX* 
2(cj2 _|_TO2a;2)3 



(4.43) 



We note that except for the small correction from the 
kinematical boundary lo = Eix oi the kj^-integral, ex- 



pression (4.43) is consistent with the Bethe-Heitler term 
derived by Zakh arov in th e dipole approximation [ [l7[ . 



Equations ( 4.42 ) or ( 4.43 ) allow to determine the only 
free parameter C from a comparison with well-tabulated 
Bethe-Heitler scattering cross sec tion [ |35|]. Once, C is 
fixed, the radiation spectrum (4.5) in the dipole approx- 
imation provides thus a parameter free prediction of the 
measured in-medium energy loss. Zakharov has pursued 
this strategy to fix C in his q^- and k^- integrated spec- 
trum ( 2.49 ), and this has lead to a very successful de- 
scription of the SLAC-146 data on radiative energy loss [ 

ITtI. 



V. THE DIPOLE PRESCRIPTION FOR QCD 

To calculate the integrated QCD radiative energy loss 
of a hard coloured parton, Zakharov [ n7| has used a 



very simple substitution in his energy loss formula (2.49). 
This QCD dipole prescription consists in replacing the 
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dipole cross section 17(0; p) in the calculation of the QED 
radiation spectrum by a combination of three dipole cross 
sections [ ^^ : 

O-QED (P, x) = cr(x p) * 

9 _ _ 1 _ 

crQCD(p,2:) = - {a{p) -V a{{l - x) p)} ~ -a{xp). (5.1) 



we have chosen a proportional to the elastic q- 

The heuristic 



cf. subsection VA 



Here 

q Mott cross section, 
argument for the prescription (5.1) starts from the repre- 
sentation of the projectile quark in the light cone frame as 
a superposition of the bare quark and higher Fock states, 



Iprojectile) = \q) + \q-f) + 



(5.2) 



If all Fock components interact with the external po- 
tentials with the same amplitude, then the coherence 
between these amplitudes is not disturbed, and no 
bremsstrahlung is generated. The radiation amplitude 
depends hence on the difference between the elastic scat- 
tering amplitudes of different fluctuations. The dipole 
cross section a{x p) contains information about this dif- 
ference since it arises from averaging in (2.2C) the (part 
of the) amplitude Aig for the "ingoing" \q) with the (part 
of the) complex conjugated amplitude Mg for the "out- 
going" higher Fock state \q^)- In the transverse plane, 
the separation of the \q"f) fluctuation from the ingoing 
q can be estimated via the uncertainty relation: x± ~f oc 



"Se cx; (1 — x)/k± and x±q (x 



ex x/Zcj 



xEiAE "^ V" ■^y/-^ ^,^^ ^j_q v^ (i-x)EiAE 

The transverse distance p is linked to the transverse mo- 
mentum in ( ^.22 ) by a Fourier transform and we may 
think of p as the transverse size of the g-7 fluctuation. 
With the transverse center of mass of this fluctuation at 
the position of the ingoing quark \q), see Figure 0, the 
transverse distance between the charged components q 
and q is then x p, and it is the dipole of this size which de- 
termines the radiation spectrum. For related arguments, 
see also Ref . [ M . 



Y 

X- 



q 



q 



(1-x)p- 



►l-^-xp--*l 



FIG. 2. Separation of dipoles in the transverse plane: the 
q-7 fluctuation of transverse size p gives rise to a QED radi- 
ation spectrum determined by the q-q dipole of size x p. In 
QCD, where the radiated gluon is charged, all three dipole 
contributions should be considered. More details are given in 
( p.l| ) and the following text. 

In QCD, the emitted gluon is charged too, and aside 
from the q-q dipole of size x p, there is a g-q dipole of 
size p and a g-q dipole of size {1 — x) p. The prefactors 
introduced for these dipoles in the prescription (5.1) stem 
from interpolating between the limits 



9 _ 

limCTQCD(p,a;) = 7 cr(p) 

x—>0 4 

liniCTQCD(p,a;) = ct(p) ■ 



(5.3) 
(5.4) 



In the limit x — + 0, the qq pair is indistinguishable from 
a pointlike color-octet charge and the ratio 9/4 of the 
octet to triplet couplings arises. In the opposite limit 
X —^ 1, the gluon-quark pair at vanishing separation is 
indistinguishable from a pointlike quark. 

The question to what extent this intuitive physical pic- 
ture leads to the correct QCD radiation cross section was 
addressed recently in Ref. [ [2^ . There it is argue d that 
Zakharov's integrated energy loss formula ( 2.4E ) coin- 
cides for QCD with the result of their calculation based 
on time-ordered perturbation theory. For a discussion of 
the terms neglected in this approach, see Ref. [|3J|. Here, 
we contribute to this discussion by comparing the trans- 
verse momentum d epe ndence resulting from the QCD 
dipole prescription (5T) to the transverse momentum de- 
pendent radiation cross sections for A^ = 1 and N — 2 
scattering centers, calculated [ M in time-ordered per- 
turbation theory. 



A. The Bertsch-Gunion spectrum for A'^ = 1 in the 
QCD dipole prescription 

In the present subsection, we test the QCD dipole pre- 
scription (5.1) for the case of A^ = 1 scattering. We start 
from the A^ = 1 QED Bethe-Heitler spectrum ( |3.14| ) in 
the form 



^5^BH-QED 

d{\nx) dq± dk^ 
-1 



Aar. 



9 ? 



' (2 7r)2 k2^(ki-.Tq. 
d(xp 



■T 



(27r) 



cTQEuix p) exp {i X q± ■ p} . (5.5) 



To specify the absolute size of cr in (5.1), we relate the 
elastic QED e^-e^ cross section to the corresponding 
QCD one by 



-1 



dr 

dr 

(2^ 



cr(r)e^'i^ 
CT(r)e*i^ 



AZ^aL 



[IVP 



-qi)2 
4 a? 



(Af2 



-,2 \2 



where Ci denotes the colour factor, 2 Ci 



1, 



(5.6) 



for 



4 

9' -^' 4' 

q-q, q-g, g-g. Also, to change to the QCD-case, we re- 
place the coupling to the emitted photon in ( ^.5| ) by 
ckom —> ctg Ca, where the Casimir Ca of the adjoint rep- 
resentation accounts for the emitted gluon. With this in- 
put, we follow the dipole prescription and substitute (5.1) 
in the Bethe-Heitler cross section (5.5). After rescaling 
qj^ -^ i qx on both sides of the equation, we find (note 
that a{p) = ct((1 - x) p) + 0{x)) 
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dVBG-QCD ^^^^ 



qi 



d(lnx) dq± dkj 



ki(k. 



■ r 



'^"^'QCD(qi) = / c?qi±dq2j 



9 "s 



(M2+q2^)2 



9 (M2 



qi)^ 



0(x), 



(5.7) 



9 "s 



This is the Bertsch-Gunion radiation spectrum [ |3g] 
times the elastic Mott cross section for a Debye-screened 
coloured potential. For A^ = 1, the QCD dipole pre- 
scription turns to leading order 0{1/E) the Bethe-Heitler 
QED radiation spectrum exactly into the Bertsch-Gunion 
QCD radiation spectrum. This success does not depend 
on the particular combination of colour dipoles intro- 
duced in (^). It is rather due to the scaling property 
of the dipole size in (pTq): replacing in the argument of 
a the transverse extension x p —f p modifies both: (i) 
the x^ rapidity distribution of the QED spectrum into 
the flat rapidity distribution of the QCD spectrum and 
(ii) the characteristic x qj_-dependence of the QED spec- 
trum into the characterisitc q^-dependence of the QCD 
radiation spectrum. 



B. Corrections to the QCD dipole prescription for 

N = 2 



The interference terms of QCD multiple scattering 
bremsstrahlung radiation are more complicated than 
their QED analogue. Especially, since the gluon is 
charged, it can rescatter on additional potentials, and 
these cascading contributions do not factorize from the 
gluon production amplitudes. This effect does not exist 
for the QED bremsstrahlung: there, the only in medium 
modification of the produced photons is the dielectric 
suppression which is clearly a final state effect which fac- 
torizes from the production amplitude. 

This qualitatively different structure of QCD rescatter- 
ing processes makes it interesting to study the accuracy 
of the QCD dipole prescription for N = 2 scatterings. 
To this aim, we compare the N — 2 radiation spectrum 
(3.15) to the result of the full N — 2 QCD radiation spec- 



trum calculated in time ordered perturbation theory in 
Ref. IM: 



,5 full QCD 
" "N=2 

d{\iix) dq_L dkj 



CaCHs 



dV? 



Bl 



Bl 



-fi^^2(12) 



Ca 



2Cf 

2 i?2 • 52(12) C0S(W2 Af) 

Bl ■ -62(12) cos(cj2o) Ai) j I , 



QCDlqi 



Bl ■ B2 cos{lui At) 



(5.8) 



0<i±-'i2±f 



LU20 — UJ2 — loq. Equa- 



whereuJo = ^,UJ2- 2<^ , 

tion (5.8) is obtained by squaring the amplitude in Eq. 



(2.20) of Ref. [ ^ which corresponds to the seven dia- 
grams in Fig. g. We have used the following shorthands: 



ft ^ h 



(M2 + q2j_)2 

k_L - qi± 



(^'^Hq-L -qi± -q2_L), 



B, 



(12) 



k^^ (k_L - qi±j" ■ 

k^ q,^ k^ qi^ q2± 



(k±-qi±)2 (kx - qi_L - q2±)^ 



(5.9) 
(5.10) 
(5.11) 



The N — 2 radiation spectrum (5.8) is obtained for the 
case of two scattering centers placed at fixed longitudinal 
positions zi — ti^ Z2 = ^2 with At = t2 — ti. In con- 
trast, the N = 2 calculation presented in section IIIB 
considers two scattering centers at arbitrary longitudinal 
positions distributed according to a homogeneous den- 
sity no within a longitudinal extension L. In fact, the 
interference terms in (5.8) depend in a different way on 
the kinematical variables than the interference terms in 
( |3.16 ). Moreover, we cannot compare the At — > limit of 
both expressions, since (5.8) is obtained in time ordered 
perturbation theory where ti < t2 is used to neglect di- 
agrammatic contributions. This leaves us for a direct 
comparison of both calculations with only the limiting 
case 



lini "^ ^^-2 

At^oo d(lna;) dq± dkj 



Cachs 



c^VQCD(qj-) 



B1+B2 + — -02(12) 



(5.12) 



which corresponds to two scattering centers at arbitrary 
large relative distance. One checks that the QCD dipole 
prescription leads for the N — 2 rescattering result ( 3.26| ) 
in the L —> 00 limit 



,5 dipole QCD 

lim J-^E^ 

L^oo d{lnx) dq± dki 



CaCIs 



dVQCB{^±) 



[Bi +^1(12) I 



2 



(5.13) 



Due to the symmetry of the qix- and q2_L-integrations, 
one can exchange Bf,^^-, -^ B'^,^^. in ( ^.13[ ). Also, for 
a direct comparison with the N — 2 case ( ^.12 ), the 



opacity factor ^ that corresponds to the probability of 
two scatterings, should be dropped. 



To discuss the differences between the result (5.13) of 
the QCD dipole prescription and the N — 2 pQCD result 
(pT2), we present in Figure H the seven Feynman ampli- 
tudes contributing to the radiation spectrum ( |5.8| ). For 
each of the diagrams, we denote by Mi the entire contri- 
bution, by Mj the part of the contribution whose phase 

f 12) 
factor depends only on ti , by M^ ' the part whose phase 

factor depends on ti and t2, etc. Based on calculations 

reported in Ref. [ E3l M , one finds 
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BI cx {Ml + M. 



(1) 



Ma 



B: 



(x{mP +M3 + mP)\ 



■^2(12) oc (Ms 



M, 



(12) 



Afr)^ 



(5.14) 
(5.15) 
(5.16) 



As can be seen from Figure 0, the terms B^ and -B| 
are the Bertsch-Gunion radiation contributions off the 
first and second scattering center respectively. The 
term i?2(i2) corresponds to gluon emission around ii and 
rescattering of the gluon at ^2 ■ Since it is the gluon and 
not the quark which scatters off the second coloured po- 



tential, the term -82(12') 



( ^.12 ). These three sets of diagrams leading to Bf, 



is enhanced by a factor ^ in 

and i?2(i2) '^^^^ ^^^'~' ^"3 seen to form the building blocks 
of an effective current in the BDMPS-approach, see e.g. 
Fig. 4of Rcf. [11. 



-^corr ^ 



Ca 
2Cf 
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(5.17) 



with which the N — 2 dipole cross section ( 5.13 ) has to 
be multiplied to regain the pQCD radiation spectrum. 
[In fact, for QCD, the kj_-integrals are divergent in the 
At — > 00 limit due to the coUinear singularity that must 
be regulated by an in-medium screening scale.] For the 
kj^- differential spectrum or for finite L, corrections are 
clearly more complicated, but the important message is 
that at least for simple few-rescattering cases, the result- 
ing expressions can be compared to pQCD results and 
the deviations can be quantified. The fact that the QCD 
dipole prescription tends to underestimate the full pQCD 
radiation cross section by a significant amount (at least 
for N = 2) may indicate that it is of use as a lower bound 
for the possible radiative energy loss. 



t ...-■■"■■ 



Mi 



t ti 



M4\-. 



t ti t2 



M; 
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ti h .-•••■ 
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(g) (^ 

tl t t2 



M. 



Mf 



t ti 



M, 



^ 



FIG. 3. The seven contributions to the TV = 2 QCD radi- 
ation spectrum in time-ordered perturbation theory. In the 
large distance limit, the Bertsch-Gunion and gluon cascading 
contributions can be identified w ith p arti cular combinations 
of these diagrams, see equations (5.14) - (5.16) and the text. 



One might have expected naively that the QCD dipole 



prescription maps the two Bethe-Heitler terms of (3.26) 
onto the two Bertsch-Gunion terms Bf and i3|. This 
is only the case for the kj_-integrated cross section (i.e., 
the energy loss) if the kinematical boundary u = x Ei is 
ignored. Then, the k^-integration allows to shift k 
k± — qi_L which changes Bf 



(12) 



Bi 



dH). 



The N = 2 low opacity expansion allows a quantitative 
check of QED-inspircd QCD calculations. The deviation 
of the limiting case (BT3) from the perturbative QCD re- 
sult is clearly non-negligible: for the k_L-integrated case, 
where the B^i^^\ and ^^-contributions have the same 
size, one can read off a correction factor 



VI. CONCLUSION 

In the present work, we have extended the KST- 
formalism [ M for quantum electrodynamics in several 
ways: (i) We have supplemented it with a regulariza- 
tion prescription which correctly removes the Weizsacker- 
Williams fields of the asymptotic states from the radia- 
tion spectrum and thus warrants the correct limiting be- 
haviour of the formalism for a set of few-scattering sce- 
narios, (ii) We have calculated the corresponding angular 
integrated energy loss and shown that it implies correc- 
tions to Zakharov's energy loss formula ( 2.49 ). (iii) We 
have derived a properly regularized radiation spectrum 
in the dipole approximation which allows for a straight- 
forward numerical calculation, as well as for further an- 
alytical studies in limiting cases. Most importantly, the 
KST-formalism describes the angular dependence of the 
radiative energy loss which was shown in section IV to 
have a very rich non-analytical structure, incorporating 
intuitive properties of in-medium propagation as e.g. the 
Moliere q_L-broadening. 

In the QCD case, the situation is much more com- 
plex due to the colour structure and the rescattering of 
the gluon. We have tested the QCD dipole prescription, 




The good 



used in various previous applications [ |15, 
against few-scattering pQCD results 
news is that this prescription recovers the A^ = 1 QCD 
Bertsch-Gunion gluon radiation spectrum. For multiple 
collisions, however, even the angular integrated distri- 
butions differ from the exact result by a sizeable factor 
which approaches 17/8 at large separation. This discrep- 
ancy is due to the increase in the effective elastic cross 
section because both the jet and the gluon undergo multi- 
ple scattering. Clearly, the QCD problem still needs fur- 
ther work and the present formalism needs to be extended 
to properly take into account non-abelian features. The 
introduction of the low opacity expansion in this paper 
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was shown to be a powerful test of any proposed exten- 
sion of QED calculations to QCD. 



ACKNOWLEDGMENTS 

We thank Boris Kopeliovich and Andreas Schafer for 
several discussions, especially about the problem of re- 
gaining the Bethe-Heitler limit from the radiation spec- 
trum of Ref. [ p^. After informing them about the 
importance of the regularization prescription discussed 



in section II B, Boris Kopeliovich noted that Alexander 



Tarasov had come to the same conclusion on the basis 
of earlier calculations of one of us. We also thank Peter 
Levai and Ivan Vitev for many discussions about radia- 
tive energy loss calculations in the BDMPS-forma lism 
and for sharing with us the unpublished result ( |5.S| ) 
whose derivation will appear in [ ^ . Finally, we thank 
U. Heinz for helpful comments regarding the manuscript. 
This work was supported by the Director, Ofhce of En- 
ergy Research, Division of Nuclear Physics of the Ofhce of 
High Energy and Nuclear Physics of the U.S. Department 
of Energy under Contract No. De-FG-02-92ER-40764. 



APPENDIX A: IN MEDIUM AVERAGES 

In this appendix, we give details of the calculation of 
in medium averages, used in the derivation of the radia- 



tion probability (2.22). We start from a set of N static 
external potentials Uq at positions (ri,Zi), 



N 



[/(x) = ^C/o(r-r„e-z.)- 



(Al) 



The in medium average (• • •) averages over the elemen- 
tary scattering centers in a volume of longitudinal exten- 
sion Az and transverse radius R according to 



"^<UU 



R^Az 



f(ri,...,rN;zi,...,ZN); 



(1> = 1 



(A2) 



We are interested in the combinations T of external po- 
tentials appearing in the squared radiation amplitude 



(2.18) 



J^[r- r'; z', z] = / cxp J ^ W,{v, r'; z' , z) I 



(A3) 



WAr,r';z',z) 



Z 

i J d^ [f/(r - r„ e - z.) - C/(r' - r„ ^ - z,)] . (A4) 



With a density n{zi,ri) of scattering centers (which is 
n = N/tt R'^ Az for the homogeneous case), the average 
i) can be rewritten in the form 



N 



exp <^ Yl ^^ 




dzi dr^n{zi,Yi)e 



w^ 



N 



(for N ~^oo), 
S= / rfz,rfr,n(z„r,) (^e^.(r.r';.'..) 



(A5) 
(A6) 



The notational shorthand S introduced here does not 
grow linearly with the number of scattering centers N. 
This is crucial for the N ^ oo limit in (A5) to work. 

So far, we have not sp ecified the remaining functional 
dependence of S in ([Aq). For homogeneous targets of 
sufficiently large transverse size, it is clear that S can only 
depend on the relative path difference p(^) = r(^) — r'(^). 
We now demonstrate for a specific model that S takes the 
form 



S[p] 



rfes(e,p(e)), 



(A7) 



where I](^, p(^)) = ■^n{£,) a[p{£_)) is defined in terms of 
the so-called dipole cross section a{p) and the longitudi- 
nal density n(^) of scattering centers, see (^-21). 

The model under consideration is the abelian version of 
the Gyulassy-Wang model which takes for the elementary 
scattering centers in (Al) Yukawa potentials with Debye 
mass M, 



C/o(x) 



47r Ixl 



(A8) 



We have included in the definition of these potentials 
an extra power in the coupling constant to take the cou- 
pling of the potential to the passing electron into account. 
In the following calculation, we assume that the average 
range 1/M of the potentials Uq is much smaller than the 
size of the medium, z' — z ^ 1/M. Also, we use that 
density fluctuations are negligible on the scale 1/M. Ex- 
panding S to leading non- vanishing order in the coupling 
constant e, we find 



a(p(0)=2 / de[((C/o(0,0)C/o(0,0)> 

-((C/o(O,0)C/o(p(O,n»], (A9) 



where 
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(([/o(0) Uo{x))) = / dx Uo{x) U„{x - x) . 



(AlO) 



The integrals in ( [A9| ) can be done analytically, leading 
to Zakharov's result f B2| 



^(P) 



M2 



■8TT{l~MpKi{Mp)) 



(All) 



In our study of the k^-dependence of the radiation prob- 
ability, we shall encounter the Fourier transform of S[p] 
with respect to the transverse momentum transfer qj^ 
from the medium to the electron, 



S(q±) 



dp 



{2n\ 



2 O-(P) exp {i q± • p] 



(M2 + q2^)2 



M2 



For n on-vanishing momentum transfer, the first term of 
( [A12| ) is the Mott cross-section for electron scattering off 
the single external potential (Al) in the relativistic limit. 
The size of the prefactor of the J^^(qj^) contribution is 
such that the q^-integral over ( A12 ) vanishes. This war- 
rants that the dipole cross section o'(p) vanishes at van- 
ishing dipole size p — Q. One can trace back the origin of 
this second term: in the expansion of (A4) to second or- 
der in the coupling constant, it arises from contributions 
which are second order in U{r — Ti) but zeroth order 
in U{r' — Vi), or vice versa. This amounts to multiply- 

(2) 

ing the second order term of Ml^ with the zeroth order 

(0)* 

term M\^ of the complex conjugated amplitude. The 

contribution Mr^' x Afi,- to the radiation cross section 
vanishes however due to energy momentum conservation. 
Calculating such a contribution from the Furry approxi- 
mation (2.5), (|2.6|), this energy momentum conservation 



appears in the form of a (5-function 6^'^'{q±), resulting in 
an additional contribution to the Mott cross section. For 
a comparison with experimental data, one has to assume 
a transverse momentum transfer q± which is always fi- 
nite but can become arbitrarily small. Then the second 
term in (A12) drops out in all physical quantities. 



APPENDIX B: EIKONAL EXPRESSIONS FROM 
PATH INTEGRALS 

In this appendix, we give details of the solutions of 
the path integrals, used for the derivation of the radi- 



ation probability {\Mfi\'^) in (2.22). We start from the 
observation that with the help of the averaging (^.20 ) 



S{z\r{z'),r'{z');z,r{z),r'{z)\p) 

EE {G{z\r{z');z,r{z)\p) G*{z' ,r'{z'); z,r'{z)\p)) 

VpVpe^py-^J 'pfj\ T[p;z',z], (Bl) 

p{z)^r{z)~r'{z), p{z) = r{z) + r'{z). (B2) 



Here, the dot denotes a derivative with respect to the 
longitudinal direction and the integral / Vp is trivial. It 
results in a (5-function which constrains all possible paths 
P(C) by the mid-point rule to a straight line 



Ps{O^P{^')^+P{^)^ 

z — z z' — z 



The final result is 

S{z',r{z'),r'{z');z,r{z)y{z)\p) 
P 



(B3) 



27ri(z'- z) 



cxp ■ 



s(^,/5.(0)rfe 



,5(2) (q^). (A12) xexp 



ip 



(r(z')-r(z))^- {r'{z')-r'{z)y 



2(z'-z) 



(B4) 



We note that this 'eikonal path' ps does not result from 
an eikonal approximation. It is rather obtained from a 
calculation on the cross section level, once the ingoing 



and outgoing states "^ p of (2^), (2^) are adopted as 



starting point: the Furry approximation is a systematic 
expansion in l/E and thus closely related to the eikonal 
approximation. 



path in M,: 




difference between 
the two patfis 



p^ patfi in M(| 



FIG. 4. Diagrammatic representation of th e path integrals 
appearing in the radiation probability (2.22). According to 



(Bl), the path differences pi and p2 at early times £, < z and 
late times ^ > 2' are fixed to straight lines. Between z and z' , 
the path difference can change due to interactions with the 
medium, but the boundary conditions are fixed. 



The explicit expression (B4) allows to simplify the ra- 



diation probability (pTTq) by using 

j dv2 dr'2 e-^P-(-^^-'-^)5(z+, r2, r^; z', p, r'jpa) 



z+ 



- exp < -ip2i ■ P2 - / S(f,P2)df| 
/"dridr;e'Pi-^(''i-'"i)5(z,r,p';z_,ri,r;|pi) 

= -exp<^ ipii -pi - / E(^,pi)rfO , 
where we have introduced the coordinates 



(B5) 



(B6) 



pi=r p 
P2 = p r' 



pi=r + p 

P2 = P + r' 



(B7) 
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The integrals (B5), (B6) depend only on the relative dis- 



We can now simplify the radiation probability {\Mfi 



tances pi, p2- Also, due to the boundary conditions, the of (2.18) with the help of (B5), (JBq) and (B15). To this 



path Ps{C} turns out to be ^-independent in (B5) and 



aim, we observe first that with a partial integration, the 
derivative d/dr d/dr' of the interaction vertex (2.17) 



Equation (pj) provides an explicit expression for the changes into —d/dpi d/dp2- Then we change to the 
in-medium average of pairs of Green's functions as long integration variables 
as both Green's functions are for particles with the same 
momentum p. To calculate the radiation probability 
( ^.17 ), we have to extend this calculation to the case of 
two different momenta. To this end, we use the identity 



T = Pl- P2 

4 4 



T = pi+ P2 



{r + p' + p + r') 



(B17) 



1 = 7V(r, p) I Vre exp < —x 



where the norm M{r, p) is given in terms of the boundary 
conditions for the path integral over r^., 



)The only r-dependence of the integrand of (|Af/ip) is in 
, (B8) the exponential of the second hne of (B15). Rewriting 
this in the new integration variables (B7) and (B17), we 
find 



2t:i {z - z) 

J\f(r,p) = ^ exp 

px 

re{z) = r{z) = J 



ipx 



2{z' - z) 



{r - pY 



(B9) 



dr exp 



ip 



(Pl - P2) • T- - f ((Pl - P2) + tY 



r,{z')^r{z')^p. (BIO) = ^![llf^^ exp 



This allows us to write 
{G{z',p;z,r\p{l~x))G*{z',r'-z,p'\p)) 
= 7V(r,p) I VrVr'VreT{r-r') 



xexpjl r [{l-x)P+xrl~r'^]\, (Bll) 

where r'(z') = r' and r'{z) — p' . Introducing the new 
coordinates 



r,(0 = (l-a;)r(C)-Hxre(0 
r,(0=r(0-re(0, 



(B12) 
(B13) 



the exponent of the path integral (BIO) reads 

{l-x)r'^ +xrl - r'2 = r^ - r'^ + (1 - x)xrl , (B14) 

and T appears with the argument J'{ra — r' + xrjY). In 
complete analogy to the calculation of (B4), this renders 
the path integral over (r^ -I- r') trivial and reduces the 
support of [ra — r') to a straight line. With the choice 
of boundary conditions (B9), we find 



{G{z\p;z,r\p2)G*{z\r';z,p'\p,)] 



2n i {z' — z) x 



X exp 



ip 



2 (z' - z) 
ifi 



Vri, exp 



(l-:r)(p-r)'-(r'-p')' 

^{^,Pf +xrb) 



ry, - 



where ^ = p[l — x) x and 

p/(e) = P2^ 



pi- 



z'-S. 



(B15) 



(B16) 



2{z' -z) 
ip{l~ x) 



(Pi " P2) 



px [_ 2(z' — z)x 

The radiation probability reads now 

(|M/,|2) = Arg I dpi dp2 db I dz I dz' 
X exp {iq{z' - z) + ipi±_ ■ px - ip2_L • Pi} 

xexpj-^ E(^,pi)rfe-^/s(C,P2)rfe 
ip (1 — a;) 



(B18) 



><r-pirip^ exp 



X / Prh exp ■ 



2(z' - z)x 
ip {\~ x)x 



{pl - P2) 



exp I- J:{£„pf+xrb) df 



(B19) 



This is proportional to a transverse area db which we 
devide out in what follows: the corresponding radiation 
spectrum is given per unit transverse area. The above 
expression can be simplified further if one observes that 



l^ {h+hYdc 



[pl - P2Y 



I rldC (B20) 



After a transformation pi -^ x pi, and p2 — > x p2, this 
allows to combine in (Bl£) the Gaussian exponent and 
the path integral into a new path integral over 



rc{0=Pm+rc{0, 
with boundary conditions 

r^z) = p/(z) = pi, 
rc{z') ^ Pf{z') =p2. 



(B21) 



(B22) 
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The resulting expression is given in (2.22) and dif- 
fers from the result of Ref. [ [2^ by the regularization 
exp{— e \z\} of the integrand only. The extension of the 
formalism to spin- and helicity dependent quantities is 



obtained, by inserting for r_pj and T^ in ( B19| ) the 
spin- and helicity dependent vertex functions ( 2.15D and 
(^.16), rather than the average (2.17). 



APPENDIX C: MULTIPLE SCATTERINGS: N=3 

For the case of A^ = 3 scatterings, the combined mo- 
mentum transfer qx from the medium is desribed by the 
folding of three elastic Mott cross sections S, 



dV^{q±) = - I dqij_ dq2± dq3± S(qi_L) S(q2±) E(q3x) 



-s^'^ im 



qi± - q2j 



q3_Lj 



(Cl) 



With the shorthands 



u„i =x (,Pi_L + qi_LJ , 

Um2 = X (pi_L + qi_L + q2±) 



u; 



ml 



'■m2 



(C2) 



2/i 2fi 

the radiation spectrum ( p. 22 ) expanded to third order 
reads 



rfV 



fi(lnx) dp± dkj 



0{al^T^) 



= Cprc5nf / dV^{q±) 



1 



2^(3) 



3! ^ ^ 
1 



1 



Um2^2 



(3) 



^ml^3 



(3) 



2^(3) 



3! 



ufZ, 



■ U2 • Um22, 



(3) 



- Ul • U„il^g 



, -7(3) , cr(3) 

+U„ii • Um2^7 + U2 • Umi^g 

+Um2 ■ UiZlf' + Ul • U2Z1I' 



r(3) 



(C3) 



The variables Zl stand again for the longitudinal inte- 
grals over phase factors. For a homogeneous medium of 
density rig and length L they read 



r(3) 



r(3) 



2Qr 



3 

m2 



■^Qm2 

r{3) _ 6 sin(L Qjni) - 6 L Q^i -f £^ Q^^ 
6 sin(L Q,n2) - 6 L Q„,2 + L^ Ql, 



r(3) 



r(3) 



3Q2Qm2 



(C4) 
(C5) 
(C6) 
(C7) 
(C8) 



(3) _ 6 sin(L Qmi) ~ 6 L Q„ii + L^ Q^i 

SQlQml 
^(3) _ L^ {siTL{LQml) - LQral) 

6Q1Q2 Qmi {Qml — Qm2) 

{sm{LQm2) - LQra2) 
(Qml — Qm2) Qm2 

„(3) _ (sin(£Qm2) - LQm2) 

m2 



Q2 (Qml — Qm2) 

(sin(LQmi) - LQmi) 

Q2 Qml {Qml — Qm2) 

,(3) _ (sin(£Q,n2) - LQm2) 

Ql {Qml — Qm2) Qm2 
{^IVLJLQml) - LQml) 
Ql Qml (Qml ^ Qm2) 

^(3) (sin(L Q.,„i) - i Qml) 



Ql Q2 Qml (Qml ~ Qm2) 

(sin(£Qm2) - LQyn2) 

2 
m2 



Ql Q2 (Qml — Qm2) Q 



(C9) 



(CIO) 



(CU) 



(C12) 



(C13) 



From this one finds again the simple limiting cases. We 
confirm the factorization limit 



dV 



L^o d(lna;) dpj^ dk 

r3 



r fixed 



0(al^T3) 



^prc Qni 



1 U2 

12 VQi Q2 

and the Bethe-Heitler limit 
d^a 



JdV'{q^), 



(C14) 



L^oo d(lna;) dpx dkx 



T fixed 



11 
12 



Ul 

Ql 



Um2 



/m2 



U2 

Q2 



Uml 
Qml 

2 



2 



Cpiegnf / dV^{q±) 



Uml 



Um2 



(C15) 



Also, we have checked in complete analogy to the case 
of A^ = 2 scatterings that the lowest order correction to 
the factorization limit (C14) is proportional to L^ times 
T^. The q^-intcgral of this correction diverges logarith- 
mically. This further corroborates the conclusions drawn 
in subsection IIIB. 
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